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Abstract. Given a sufficiently large universe, numerous civilizations almost surely exist. 
Some of these civilizations will be short-lived and die out relatively early in their develop- 
ment, i.e., before having the chance to spread to other planets. Others will be long-lived, 
potentially colonizing their galaxy and becoming enormous in size. What fraction of civ- 
ilizations in the universe are long-lived? The "universal doomsday" argument states that 
long-lived civilizations must be rare because if they were not, we should find ourselves living 
in one. Furthermore, because long-lived civilizations are rare, our civilization's prospects for 
long-term survival are poor. Here, we develop the formalism required for universal doomsday 
calculations and show that while the argument has some force, our future is not as gloomy 
as the traditional doomsday argument would suggest, at least when the number of early 
existential threats is small. 
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1 Introduction 

The Doomsday Argument [1, 2] 1 traditionally runs as follows. Perhaps our civilization will 
soon succumb to some existential threat (nuclear war, asteroid impact, pandemic, etc.) so 
that the number of humans ever to exist is not much more than the number who have 
existed so far. We will call such a civilization short-lived and the total number of humans 
in it Ns- Alternatively, we might survive all such threats and become long-lived, potentially 
colonizing other planets and eventually generating a large total number of individuals, Nl- 
For simplicity we will consider only two possible sizes. We expect Nl S> Ns- The ratio 
R = Nl/Ns could easily be as large as a billion. 

We don't know our chances of being short- or long-lived, but we can assign some prior 
belief or credence P(S) and P(L) = 1 — P(S) in these two possibilities. These confidence 
levels should be based on our analysis of specific threats that we have considered and the 
possibility of other threats of which we are not yet aware. 

Suppose that you hold such confidence levels at a time when you don't know your own 
position in the entire human race. Now you discover that you are one of the first Ns humans 
to be born. We will call this datum D. If the human race is to be short-lived, D is certain. 
If the human race is to be long-lived, assuming that you can consider yourself a randomly 
chosen human [5-8], the chance that you would be in the first Ns is only 1/R. Thus you 
should update your probabilities using Bayes' Rule, to get 

P(L\D) = P{L),R = P{L) < (1 1) 

11 ; P(L)/R + P(S) P(L)+P(S)R P(S)R' { ' } 

Since it is clear that we do face existential threats, P(S) is not infinitesimal. Thus P(S)R 3> 
1, P{L\D) -C 1, and doom (our civilization ending soon rather than growing to large size) is 
nearly certain. 



1 Gott [3] and Nielsen [4] make similar arguments. 
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Many counterarguments have been offered (e.g., [7, 9-12]). The specific issue which 
will concern us here is the possibility that our universe might contain many civilizations. In 
that case, we should consider ourselves to be randomly chosen from all individuals in that 
universe or multiverse. Before taking into account D, our chance to be in any given long- 
lived civilization is then higher than our chance to be in any given short-lived civilization by 
factor R. Taking into account D simply cancels this factor, so the chance that we are in a 
long-lived civilization is just the fraction of civilizations that are long-lived. (More compactly, 
since each civilization contains N$ individuals who observe D, this observation provides no 
reason to prefer one type of civilization over another.) 

Thus if there are many civilizations, the doomsday argument is defeated. However, it 
returns in another form, called the universal doomsday argument [13, 14]. We are more likely 
to observe D in a universe in which most civilizations are short-lived. While we can no longer 
conclude anything unique to our own civilization, we can conclude that most civilizations are 
likely to be short-lived, and thus ours is likely to be short-lived also. 

This paper analyzes the universal doomsday argument. In contrast to the traditional 
doomsday argument, which in almost all circumstances makes doom nearly certain, for many 
reasonable priors the universal argument gives only a mildly pessimistic conclusion. However, 
for other priors, the conclusion of the universal argument can be quite strong. 

The analysis of many civilizations in the universe can be extended to analyze possible 
civilizations that might exist according to different theories of the universe. If we take 
ourselves to be randomly chosen among all observers that might exist [9, 12], the doomsday 
argument is completely canceled, and after taking into account D we find no change in our 
prior probabilities for different lifetimes of our civilization. This assumption, equivalent to 
the self-indication assumption [9], is a controversial one, and the authors of the present paper 
are not in agreement about it. However, for the purpose of the present work we will consider 
the consequences of denying this idea, and consider ourselves to be randomly chosen only 
among those individuals who exist in our actual universe. 

In the next section we set up the formalism for calculating the probability P(L) for a 
civilization to be long-lived. To simplify the discussion, we focus on a special case where 
civilizations can have only two possible sizes; the general case is discussed in the Appendix. 
In Sections 3-6 this formalism is applied to find P(L) for several choices of the prior. Our 
conclusions are summarized and discussed in Section 7. 

2 The fraction of long-lived civilizations in the universe 

Assume the universe is large enough such that numerous civilizations exist. Assume further- 
more that a fraction fa of the civilizations are long-lived, containing Nl individuals each, 
and the remaining fraction 1 — fa of the civilizations are short-lived, containing only Ng 
individuals each. (The general case in which civilizations may have any size, rather than just 
the specific sizes Ng and Nl, is discussed in the appendix.) We will take the universe to be 
finite or the problems that arise in infinite universes to have been solved, so that the fraction 
fa is well defined. 

We do not know what fa is, brit. we know it must be in the interval, [0, 1]. We can 
represent our prior belief, or "best guess" for different values of fa, with a density function 
P(fa)) so that P(fL)dfL is our prior probability for the fraction of long-lived civilizations to 
lie within an infinitesimal interval dfL around fa. We only consider normalized priors so that 

./;,'/'(.//)'//•/. i- 
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Figure 1. The likelihood function, P(D\f), given R = 10 6 , R = 10 9 , and R = 10 



Whatever we take as our prior, we should update it after observing new evidence or 
data. The rest of this paper is concerned with updating P(/l) after considering the datum 
D: we were born into a civilization that has not yet reached long-lived status. Bayes' Rule 
gives 

P(MD)= P WM P M (2 1) 

nfLlIJ) fdf L P(D\f L )P(f L )- [2A) 

With a given f^, the probability of observing D, P(D\fi,), is just the ratio of the number 
of observers making observation D to the total number of observers. Factoring out the total 
number of civilizations, we divide the number of observers in each civilization finding D, 
which is just Ns, by the average number of observers per civilization, N^fi + N$(l — f£)- 
Thus 

In Fig. 1, we plot the likelihood function given R = 10 6 , R = 10 9 , and R = 10 12 . 
Plugging the likelihood into (2.1), 

P(h\D)= /^>/( 1 + ^ fl - 1 » . (2.3) 

The posterior distribution, P{fi\D), expresses our updated belief in different values of fi 
after considering that our civilization is not yet long-lived. 

Notice that P(fz\D) oc P(D\fi,)P(fi,) and that P{D\fi) is much larger for low values 
of Jl than for high values of fi (see Fig. 1). Therefore, when updating our prior, we place 
more credence in low values of fi and less credence in high values of Jl , and so become more 
pessimistic about the fraction of civilizations that reach long-lived status. 

We can now compute the probability that our civilization will eventually be long-lived. 
For any given fi, this probability is just fi, the fraction of civilizations that are long-lived. 
Without considering D, we would just take an average of the possible Jl weighted by the 
prior, 

P(L)= [ df L f L P(f L ). (2.4) 
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To get the posterior chance, we integrate over our posterior probability distribution for f^, 



J Jo dfLP(fL)/(l + fh{R ~ 1)) 

Given R and the prior distribution P(f£), (2.5) allows us to determine our civilization's 
prospects for long-term survival. 

Equation (2.5) is exact, but since we are interested only in R S> 1, we can always 
approximate 

p(l\d) ~ & dfL + fhR)) - dfL fLPihyjR- 1 + h) (26) 

tidf L P(f L )/(l + f L R) fidf L P(f L )/(R-i + f L ) 

In many cases the contribution from fi < to the numerator is not significant, and we 
can thus approximate the numerator by dfiPifL) = 1, giving 



P(L\D) 



df L ~ 



o 



(2.7) 



In some cases there is a cutoff on that keeps it above and in such cases we can write 



P(L\D) 



1 at 

«/l — - t — 

JL 



-1 



(2-« 



From (2.8) we can understand the general effect of the universal doomsday argument. 
The chance that our civilization will survive to large sizes is small when the integral is large. 
This happens whenever there is a possibility of small fi where the prior probability of those 
/l is large compared to the /l themselves. So, for example, if our prior gives a collective 
probability of 10 -6 to a range of fi near 10 -9 , i.e., we think there's one chance in a million 
that only a billionth of all civilizations grow large, then our chance of survival is no more 
than f0~ 3 . 

The effect of using (2.7) instead of (2.8) is that the above argument does not apply to 
fl below R . The doomsday argument is able to increase a prior probability by factor R 
at most, so scenarios with prior probabilities below R^ 1 will never be important. 

In the rest of this paper, we will consider several specific priors and see the conclusions 
to which they lead. 

3 The uniform prior 

The results for P{L\D) will depend on our prior P(/^). Let us start by taking the simplest 
prior, 

P(f L ) = l, (3.1) 

so that the fi is equally likely to have any value. This appears to be a reasonable choice of 
prior when we do not have much quantitative information on the existential threats that we 
are facing. The posterior density is then 
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Figure 2. Cumulative priors and posteriors with R = 10 6 and equal prior credence given to all ft, 

In Fig. 2, we plot the cumulative of the prior and posterior, P(/l < /) and P(/l < f\D), 
given R = 10 6 . Although the probability of low values of /l increases after considering D, 
the result is not extreme. For example, the probability that Jl > 1/2 is 5% — a low but 
non-negligible amount. 

The probability that our civilization is long-lived, after considering D, is 



P(L\D) 



1 



1 



In R R-l 



(3.3) 



Because P(L\D) ~ 1/lni?, our long-term prospects are not too bad. For example, when 
R = 1 million, our civilization's chance of long-term survival is approximately 7%. 

We can compare these results with the traditional doomsday argument. Using the 
uniform prior, the prior chance that our civilization would be long-lived is 1/2, and the 
posterior chance about 1/lni?. If we took a prior chance of survival P(L) = 1/2 in the 
traditional doomsday argument, (1.1) would give our chance of survival as only 1/{R+ 1). 
Thus, at least in this case, taking account of the existence of multiple civilizations yields a 
much more optimistic conclusion. 

We can reproduce the traditional doomsday argument even in the universal setting, 
merely by asserting that all civilizations have the same fate, so the benefit of multiple civi- 
lizations is eliminated. This would imply giving no credence to any possibilities except fi = 
(all civilizations short-lived) and /l = 1 (all civilizations long-lived). Thus we could write 



PUl) = P(L)S(f L ) + (1 - P(L))5(f L - 1) 



(3.4) 



Using (3.4) in (2.5) reproduces (1.1). 



4 N existential threats 

Of course we know that civilizations face more than one existential threat. So let us consider 
the case where there are N statistically independent threats and take a uniform prior for the 
fraction of civilizations that survive each one. Denote the fraction of civilizations surviving 
the i-th threat fi . The fraction of civilizations that survive all threats is then 



fi — /1/2/3 • • • f. 



N ■ 



(4.1) 
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and our prior for each threat is 

P(fi) = l, (4.2) 

so that 

P(h, h, ■ ■ ■ , In) = P{fi)P{h) ■ ■ ■ PUn) = 1 • (4.3) 
We can determine the density function for the overall prior, P(/l), as follows. Let I = 
| In /l I = — In fi = I i where Z j = | In I j | . Then 

P(l) = PUl)^ = P(fL)e~ l . (4.4) 

Similarly P(h) = e~ li , so I is the sum of N independent and exponentially distributed random 
variables, and P(/l) is thus given by an Erlang (Gamma) distribution, 

l N - l e~ l 

p(,) = (iv^ < 4 ' 5) 

giving 

p ^=W- (46) 

The cumulative of the prior, P(fi < /), is shown in Fig. 3(a) for N = 1 to N = 5. Note 
that our civilization's prospects for long-term survival become rather bleak as ./V increases. 
Even without considering our datum, P{L) = 1/2 N . 

Now, considering our datum D, the posterior density is, 

P(fL\D)= (1 - R) |ln/ ^ |jV " 1 (4 7) 

ULl 1 U N (1 -R)(N- 1)!(1 + f L (R - 1)) ' { ■ 1 

where Li^r is the polylogarithm function, given by 

1 f 00 t^" 1 
T(N)J z^J-l 

The cumulative of the posterior, P{fi < f\D)i is shown in Fig. 3(b) for N = 1 to N = 5 and 
R = 10 6 . 

In Fig. 4, P(L\D) is shown as a function of R for N = 1 through N = 5. After 
considering D, our civilization's prospects for long-term survival are 



U N (z) = —l __ lt - . (4.8) 



When lni? 3> JV, we can approximate 

AT! 
(ln#) J 

Notice that when N = 1, the prior distribution for ji, is uniform and the results are 
the same as in the previous section. As seen in Fig. 4, increasing the number of existential 
threats, N, decreases the probability that our civilization will be long-lived. The chance of 
survival before taking into account D is P(L) = 2~ N . Updating adds the additional factor 

P(L) {lnR) N ' 1 ; 

which is small whenever lni? S> N. However, only a power of lni?, rather than R itself, 
appears in the denominator, so the effect is more benign than in the traditional doomsday 
argument of (1.1), so long as is fairly small. 
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Figure 3. The cumulative of the (a) priors and (b) posteriors of /l given N existential threats and 
i?=10 6 . 




Figure 4. The probability that our civilization will be long-lived as a function of R given N existential 
threats. 

5 An unknown number of threats: Gaussian distribution 

It would be foolish to imagine that we know of all existential threats. For example, before the 
1930's no one could have imagined the threat of nuclear war. So there is some uncertainty 
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about the number of threats, N. To make a simple model of this effect, let us assume that a 
fixed fraction q of civilizations survive each threat, and thus fx = q N survive them all. We 
will take a Gaussian prior for the distribution of N, 



P N = \e 



-(N/N f 



where 



2 = E 



N=0 



-(N/N ) 2 _ + 1 

2 



(5.1) 



(5.2) 



For the purposes of the present section it is a good enough approximation to drop the 1 and 
just use Z = ^/ttNq/2. 

The prior expectation value for TV is 



(N) = Y^NP N ^ / NP N = -±. 

Jo V 71 " 

The prior chance that our civilization will grow large is 



(5.3) 



(5.4) 



providing that AT 1 lnq\ » 1, so that the Gaussian does not decline until q N is already small. 
Thus the prior chance of survival is quite appreciable. 

The posterior chance of survival may be approximated 



P{L\D) 



Pn 

^ R- 1 + f. 



= z 



^ R~ l +q 



N 



(5.5) 



As N increases, the numerator in the sum decreases more and more quickly. We can approx- 
imate that the denominator decreases by a factor q, until N = In lnq\, at which point 
it becomes constant. Thus for N < lni?/| lnq\, the iVth term in the sum is larger than the 
previous term by factor of order 



e -2N/N^ q = e \lnq\-2N/N^ 



(5.6) 



so the sum is dominated by terms near N = Nq \ lng|/2 or near N = lni?/| lnq\, whichever 
is smaller. 

In the former case, it is interesting to note that N is proportional to Nq, not Nq as 
one might have thought. The doomsday argument partly cancels the Gaussian suppression 
of the prior probability, which drops rapidly when N > Nq. 

We can approximate the sum by an integral, 



f 

Jo 



" dN exp ( (^-^o 2 |lng|/2) 2 + N§\ Inqf 



so 



^o 2 



P(L\D)^ - exp 



4 



^/ttNq exp 



1 -0.12JVg 

2 



N^\lnq\ 



(5.7) 
(5.8) 
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Figure 5. The probability that our civilization will be long-lived as a function of R with the Gaussian 
prior of (5.1) for the number of threats, with the parameter iVp = 2, 4, 6, 8, 10. 

where the last step is for q = 1/2. If iVo is small, this is not too pessimistic a conclusion. 
For example with N = 5 threats we find P{L\D) « 0.02, but for N = 10 we find P(L\D) « 
2 x 10~ 6 . 

To be in this regime, we require that Nq\ ]nq\/2 < lni?/| lng|, i.e., 



for q = 1/2. For iVo = 10 this requires R > 2 x 10 10 , which is not unreasonable. 

Fig. 5 shows our prospects for long-term survival as a function of R for various values of 
No. Because (N) « No/\/tt rs No/2, a given No corresponds roughly to choosing N = No/2 
in the previous section. 

Comparing Fig. 5 for each iVo to Fig. 4 for N = 2No, we see that, as R becomes larger, 
Fig. 4 shows lower survival chances than Fig. 5. In the present case the doomsday argument 
amplifies the probability of small / that arise from a large number of threats, but because 
we took fixed q we don't consider threats that have very low survival probability. In the 
previous section, we fixed N, but the uniform prior gave some chance to arbitrarily small /. 

6 An unknown number of threats: exponential distribution 

The reason that the conclusion of the previous section was not too pessimistic is that the 
Gaussian prior for P/v very strongly suppressed numbers of threats N S> iVo- If we take a 
distribution that falls only exponentially, the result may be quite different. So consider now 



In 12 > 



Ni\ In q 



0.24iV n 2 



(5.9) 



2 



P N = (1 - s)s 



N 



(6.1) 



for some s < 1. The prior expectation value for N is 




(6.2) 



and the prior chance of survival is 




(6.3) 
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The posterior chance of survival is given by 



P(L\D) 



^ s N 



^ R-i+q 



N 



-1 



(6.4) 



If s < q, our prior credence in N threats decreases faster than the chance of surviving TV 
threats, and the terms in the sum in (6.4) decrease. Except when q is very close to s, we can 
ignore R and get an optimistic conclusion 

P(L \ D) ^ . (6 .5) 

But if s > q, the situation is different. In this case, our credence in a large number of 
threats is higher than the survival chance, and then the universal doomsday argument acts 
to increase our prior for facing many threats and so decrease our expectation of survival. 
The terms in the sum in (6.4) increase as (s/q) N until we reach some N where the existence 
of R begins to matter. This happens when q N ~ -R -1 , i.e., when 

N « N' = ^ . (6.6) 
I l n( ?l 

For N > N' we can ignore q to get Rs N in the sum of (6.4). Thus we can split our sum into 
two parts, which give 

s/q-1 + l-s !-V ' 

and so 

P(L\D)f* ^Z±R\^s\/\lnq\-l _ ( g_g) 

s — sq 

As long as s is significantly larger than q, this is very small. As an example, we can choose 
q = 1/2 and s = 3/4. This gives the prior probability P(L) = 0.8 here, while (5.4) and (5.2) 
give P(L) = 0.81. While these are nearly the same, (6.8) with R = 10 9 gives 

P{L\D) ss 4 x 10" 6 , (6.9) 

must lower than P(L\D) « 0.02 from (5.8). 

Fig. 6 shows our prospects for long-term survival as a function of R for various values 

of s. 



7 Summary and discussion 

Bayes' Theorem tells us how the probabilities we assign to various hypotheses should be 
updated when we learn new information. The Doomsday Argument is concerned with the 
impact of the important piece of information that we called D — that we are among the 
first Ns humans to be born. Earlier investigations [1-4, 13, 14] suggested that the resulting 
probability for our civilization to be long-lived is suppressed by a huge factor R = Nl/N$ S> 
1, where is the size a civilization may reach if it does not succumb to early existential 
threats. Here, we attempted a more careful analysis by considering a number of possible 
choices of prior probabilities. We found that, with a seemingly reasonable choice of the prior, 
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Figure 6. The probability that our civilization will be long-lived as a function of R with the ex- 
ponential prior of (6.1) for the number of threats, with the parameter s — 1/2,2/3,3/4,4/5,5/6 
corresponding to (N) being 1, 2, 3, 4, 5. 

our chances of long-term survival are suppressed by a power of lni?, rather than by R, with 
the power determined by the number of threats N. If N is not too large, the probability of 
long-term survival is about a few percent. 

This conclusion has been reached by assuming a flat prior distribution, P(fi) = 1, 
for the fraction of civilizations fi surviving statistically independent threats (labeled by i). 
This appears to be a reasonable assumption, reflecting our present state of ignorance on the 
subject. We also considered a prior where the survival probability for each threat is a fixed 
number q, while the number of threats N is a stochastic variable with a Gaussian distribution 
of width Nq. Once again, we find, assuming that q is not too small and Nq is not too large, 
that the Bayesian suppression factor is much smaller than suggested by the naive Doomsday 
Argument. 

In our analysis, we adopted a model where civilizations can have only two possible 
sizes, Ns and Njj. This model is of course not realistic in detail, but it may well capture 
the bimodal character of the realistic size distribution. Civilizations that actively engage 
in colonizing other planetary systems and reach a certain critical size are likely to grow 
extremely large, while civilizations confined to their home planet must be rather limited in 
size. 

Even though we found a greater survival probability than in Refs. [1-4, 13, 14], our 
conclusions can hardly be called optimistic. With the priors that we considered, the fraction 
of civilizations that last long enough to become large is not likely to exceed a few percent. 
If there is a message here for our own civilization, it is that it would be wise to devote 
considerable resources (i) for developing methods of diverting known existential threats and 
(ii) for space exploration and colonization. Civilizations that adopt this policy are more likely 
to be among the lucky few that beat the odds. Somewhat encouragingly, our results indicate 
that the odds are not as overwhelmingly low as suggested by earlier work. 

A Arbitrary possibilities for civilization size 

In this appendix we give a formalism for considering all possible civilization sizes. Let a 
scenario be a set of numbers f n , n = . . . oo giving the fraction of civilizations that have 



- 11 - 



each size n. Since every civilization has some size, we must have Yl n /n = 1- We will write 
f for the entire vector of numbers f n . The average number of observers per civilization in 
scenario f is 

n(f)=J>/„, (A.l) 

n 

which we will assume is finite. 

Now let P (f) denote the prior probability that we assign to each possible scenario f. 
The probabilities must be normalized, 

JdSP(f) = l, (A.2) 

where 

f (if denotes f df Q [ dfx f df 2 . . ■ . (A.3) 
J Jo Jo Jo 

We suppose that P(f ) already contains a term such as (5(1 — ^ n f n ) that excludes unnormal- 

ized f . 

Let no denote the size of our civilization. We will not be concerned here with issues 
related to civilizations with fewer than uq observers, so let us suppose that P(f) is supported 
only when f n = for n < hq. Now we consider the datum D, that we are in the first no 
individuals in our civilization. The chance for a randomly chosen observer to observe D is 

< A - 4 > 

Now let A be some property of f , such as the average size of a civilization or the chance that 
the civilization has more than a certain number of members. The average value of A not 
taking into account D is 

(A) = [ d£A(f)P(f). (A.5) 



Now we take D into account using Bayes' Rule. We find 

P(t\D) = P(D|f)P(f) = (A 6) 

U ; f df' P(D\f')P(P) J df' P(f')/n(f')' 1 

which is the arbitrary-size generalization of (2.3). The average value of A taking into account 
D is 

J dfA(f)P(f)/n(f ) 
/dfP(f)/n(f) 

A particularly simple case is when A({) = n(f). The expected value of the total size of 
our civilization taking into account D is 



^ D= rwf pm/ m - (A - 8) 
J diP{f)/n{i) 

Alternatively, let /z,(f) = Ylri=N L fn ^ e ^ ne f rac ti° n of civilizations in scenario f that 
grow larger than some threshold N^. The posterior chance that our civilization will reach 
this threshold is then 

P(L\D) -(fr)\D- fXfLWPWWf) (A g) 

P(L\D) - {jl)\D - /dfp(f)/n(f) • (A.9) 
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This is just the fraction of large civilizations in the different scenarios, weighted by prior 
probability of the scenario and the inverse of the average civilization size according to that 
scenario. 

Unfortunately, the set of possible priors is so large here that it is difficult to make any 
progress. It seems likely to us that civilizations will either remain confined to a single planet 
and eventually be wiped out by some disaster, or spread through the galaxy and grow to 
large size. We can approximate this by considering P(f) supported only at two sizes N$ and 
N L , 

P(f) = P(f L )S(f s + f L -l) J] 5(f n ), (A.10) 

n^N L ,N s 

where we have written fi for f^ L and fs for /at s . Then we recover the results in the main 
text. Integrals P(f)df become P(f L )df L , n(f) is just (1 - f L )N s + N L f L , f L (f) is f L , (A.6) 
becomes (2.3), and (A. 9) becomes (2.5). 
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